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ABSTRACT 

We explore the relationship between the metallicity of the intracluster gas in clusters 
of galaxies, determined by X-ray spectroscopy, and the presence of cooling flows. Us- 
ing ASCA spectra and ROSAT images, we demonstrate a clear segregation between 
the metallicities of clusters with and without cooling flows. On average, cooling-flow 
clusters have an emission- weighted metallicity a factor ~ 1.8 times higher than that 
of non-cooling flow systems. We suggest this to be due to the presence of metallicity 
gradients in the cooling flow clusters, coupled with the sharply peaked X-ray surface 
brightness profiles of these systems. Non-cooling flow clusters have much flatter X- 
ray surface brightness distributions and are thought to have undergone recent merger 
events which may have mixed the central high-metallicity gas with the surrounding 
less metal-rich material. We find no evidence for evolution in the emission-weighted 
metallicities of clusters within z ~ 0.3. 

Key words: galaxies: clusters: general - cooling flows - intergalactic medium - 
X-rays: galaxies 
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1 INTRODUCTION 

X-ray spectroscopy provides an accurate measure of the 
metallicity of the hot intracluster medium (ICM) in clus- 
ters of galaxies. The strengths of the various emission lines 
relative to the continuum reveal the abundances of the emit- 
ting elements relative to hydrogen. Such measurements are 
particularly clear for the case of iron in rich clusters where 
the K-shell lines are typically well-defined in the X-ray spec- 
tra. The mass of the ICM dominates over the visible mass 
in stars in a cluster by a factor of 2-5 (Arnaud et al. 1992; 
David, Jones & Forman 1995; White & Fabian 1995). Metal- 
licity measurements from X-ray observations thus provide 
firm constraints on the history of metal production within 
the cluster potential wells. 

For low-redshift clusters, the observed abundance of 
iron in the ICM is approximately 1/3 solar {e.g. Edge & 
Stewart 1991). The strong correlation between the mass of 
iron and the total optical light from elliptical and lenticular 
galaxies within a cluster suggests that early-type galaxies are 
responsible for the bulk of the enrichment (Arnaud 1992). 
Mushotzky et al. (1996) showed that the relative abundances 
of individual elements such as Si, S and Fe determined from 
X-ray spectra suggest that most of the metals originate from 
type II supernovae. Mushotzky and Lowenstein (1997) fur- 
ther demonstrated that the iron abundance in rich clusters 
shows little evolution between z 0.3 and now, suggest- 
ing that most of the enrichment of the ICM occurred at 



high redshifts {z > 0.3). This is consistent with current 
semi-analytic models of galaxy formation (e.g. Kauffmann 
& Chariot 1997) which find that more than 80 per cent of 
the metal enrichment occurs at z > 1. 

X-ray observations of clusters of galaxies show that in 
the central regions of most (70 — 90 per cent) clusters the 
cooling time of the ICM is significantly less than the Hubble 
time (Edge et al. 1992; Peres et al. 1997). The observed 
cooling leads to a slow net inflow of material towards the 
cluster centre; a process known as a cooling flow (Fabian 
1994). The X-ray imaging data show that gas typically 'cools 
out' throughout the central few tens to hundreds of kpc in 
the clusters. Recent spatially resolved X-ray spectroscopy 
has confirmed the presence of distributed cool (and rapidly 
cooling) gas in cooling fiows, with a spatial distribution and 
luminosity in excellent agreement with the predictions from 
the imaging data (Allen & Fabian 1997). The cooling flow 
in an X-ray luminous cluster can account for up to --^ 70 per 
cent of the total bolometric luminosity of the system (about 
half of this luminosity being due to material cooling out of 
the flow, and the rest due to the gravitational work done on 
the gas as it flows inwards; Allen et al. 1998). 

Where abundance measurements have been compiled 
for large samples of clusters, a significant dispersion in the 
metallicity (which is primarily determined from the iron 
abundance) for clusters of a fixed X-ray luminosity and/or 
temperature has been revealed. In particular, the metallicity 
appears to depend on whether or not a cluster has a cool- 
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ing flow (Yamashita 1992; Fabian et al. 1994a) in the sense 
that cooling-flow clusters tend to have higher metallicities. 
Fabian et al. (1994a) speculated that this may be due to dif- 
ferences in the origins of the clusters, to the intracluster gas 
in the cooling-flow clusters being inhomogcncous with the 
cooler gas being more metal-rich, or due to the presence of 
metallicity gradients in these systems. Those clusters with- 
out cooling flows are assumed to have experienced a major 
merger event which has mixed the gas, loading to a lower 
emission-weighted metallicity. 

In this paper we present metallicity measurements for a 
sample of 30 X-ray luminous clusters observed with ASCA 
and ROSAT (sec Allen et al. 1998 for details). Wc identify 
the cooling flows in the sample and explicitly account for the 
effects of the cooling gas on the X-ray spectra. We argue that 
an abundance gradient in the cooling flow clusters under- 
lies the variation seen in the emission-weighted metallicities. 
Throughout this paper, we assume i?o=50 kms~^ Mpc~^, 
f2 = 1 and A = 0. 



2 OBSERVATIONS AND DATA ANALYSIS 

Full details of the reduction and analysis of the total sam- 
ple of 30 clusters are given in Allen et al (1998). ROSAT 
High Resolution Imager (HRI) images were used to map the 
X-ray surface brightness profiles of the cluster cores and de- 
termine whether the individual clusters have cooling flows 
or not. The ASCA spectra were then fitted with a series 
of appropriate spectral models (see below) from which the 
metallicity measurements wore made. 

For the purposes of this paper, we have classified the 
clusters into subsamples of cooling-fiow (CF) and non- 
cooling flow (NCF) systems. CFs are those clusters for which 
the upper (90 per cent confidence) limit to the central cool- 
ing time, as determined from a dcprojcction analysis of the 
ROSAT HRI X-ray images, is less than 10^° yr. NCFs are 
those systems with upper limits to their central cooling times 
> 10^° yr. [The 'central' cooling time is the mean cooling 
time of the cluster gas in the innermost bin included in the 
dcprojcction analysis, which is of variable size (sec Allen et 
al. 1998 for details). The use of a fixed physical size of 100 
kpc for the central bin leads to very similar results (see Ta^ 
ble).] Using this simple classification we identify 21 CFs and 
9 NCFs in our sample. The mean redshift for the subsamples 
of both CF and NCF clusters is z = 0.21. 

The modelling of the ASCA spectra was carried out us- 
ing the XSPEC spectral fitting package (version 9.0; Arnaud 
1996). The spectra were modelled using the plasma codes 
of Kaastra & Mewe (1993; incorporating the Fe L calcula- 
tions by Liedhal in XSPEC version 9.0) and the photoelec- 
tric absorption models of Balucinska-Church & McCammon 
(1992). The data from all four ASCA detectors were anal- 
ysed simultaneously with the fit parameters linked to take 
the same values across the data sets. The exceptions were 
the emission measures of the ambient cluster gas in the four 
detectors which, due to the different extraction radii used, 
were allowed to fit independently. 

Three models were fitted to the spectra. Model A, con- 
sisted of an isothermal plasma in coUisional equilibrium, at 
the optically-determined redshift for the cluster, and ab- 
sorbed by the nominal Galactic column density (Dickey & 
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Figure 1. The emission- weighted metallicity of the cluster gas 

(determined with spectral model A) versus the central cooling 
time (from the dcprojcction analysis of the ROSAT HRI images). 
CFs are plotted as filled circles and non-cooling flows as open 
symbols. The figure illustrates a clear segregation in the mean 
emission-weighted metallicities for CF and NCF clusters. 

Lockman 1990). The free parameters in this model were the 
temperature (kT) and metallicity (Z) of the plasma and 
the emission measures in the four detectors. (The metallic- 
ities axe determined relative to the solar values of Anders 
& Grevesse (1989) with the different elements assumed to 
be present in solar ratios.) Secondly, model B, which was 
identical to model A but with the absorbing column density 
(A'h) also included as a free parameter in the fits. Thirdly, 
model C, which included an additional component explic- 
itly accounting for the emission from the cooling flows in 
the clusters. The material in the cooling flows is assumed to 
cool at constant pressure from the ambient cluster temper- 
ature, following the prescription of Johnstone et al. (1992). 
The normalization of the cooling-flow component is param- 
eterized in terms of a mass deposition rate, M, which was 
free parameter in the fits. (The metallicity of the cooling 
gas was assumed to be equal to that of the ambient ICM.) 
The cooling flows were also assumed to be absorbed by an 
intrinsic column density, AiVn (Allen & Fabian 1997 and ref- 
erences therein), which was a further free parameter in the 
fits. The metallicity of X-ray absorbing material was fixed 
at the solar value (Anders & Grevesse 1989). 



3 THE METALLICITY MEASUREMENTS 

The metallicity measurements for the clusters axe summa- 
rized in the Table. For the CF clusters we list the results 
obtained both with spectral models A and C {i.e. with and 
without the emission from the cooling fiows accounted for 
in the spectral analysis). For the NCFs only the results for 
model A axe listed. 
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Table 1. Summary of the cooling time, metallicity and baryon fraction measurements 
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NOTES: Column 2 lists the redshifts for the clusters. Columns 3 and 4 list the central cooling times and the mean cooling times within 
100 kpc of the cluster centres determined from the deprojection analysis of the ROSAT HRI images. Columns 5 and 6 summarize the 
metallicity measurements for the clusters from the ASCA data using both the isothermal spectral model (A) and, for the CF clusters, 
the more sophisticated model incorporating the cooling- flow emission component (model C). Column 7 lists the ratios of the X-ray gas 
mass to the total mass at a radius of 500 kpc in the clusters. The errors on the central cooling times are the 10 and 90 percentile values 
determined from 100 Monte Carlo simulations in the deprojection analysis. Errors on the metallicity results and baryon fractions are 90 
per cent {Ax^ = 2.71) confidence limits. Errors on the mean values are the standard deviations of the distributions. 
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The results listed in the Table demonstrate a clear seg- 
regation between the metallicities of the CF and NCF clus- 
ters. For the CF systems, the mean metallicity determined 
with spectral model A is 0.34. For the NCFs, the mean value 
determined with the same spectral model is 0.21. The ap- 
plication of a Students t-test (accounting for the possibiUty 
of unequal variances in the two distributions; Press et al. 
1992) indicates that the mean metallicities for the CF and 
NCF clusters (determined using spectral model A) differ at 
^ 99.9 per cent confidence. (The probability that the means 
of the two distributions are equal is 1.1 x 10~^.) Note, how- 
ever, that an F-test (Press et al. 1992) indicates the vari- 
ances of the two distributions to be consistent at the 91 per 
cent confidence level. 

For the CF clusters, the metallicities determined with 

spectral model C, which incorporates the cooling-flow com- 
ponent, are slightly higher {Zcf = 0.38) than those inferred 
with spectral model A. The use of the cooling-flow model, 
where appropriate, therefore only enhances the discrepancy 
between the metallicities of the CF and NCF clusters. (The 
probability that the moan metallicities for the CF clusters, 
determined with spectral model C, and for the NCFs, deter- 
mined with spectral model A, are equal is only 4.3 x 10~^). 
We note that the use of spectral model B in place of model 
A (i. e. allowing the absorbing column density acting on the 
clusters to be a free parameter in the fits) leads to very 
similar results. 

The mean metallicity for the whole sample of 30 clus- 
ters, determined with spectral model A, is Z = 0.303±0.085. 
If spectral model C is used for the CF clusters, the over all 
sample mean rises to Z = 0.329 ±0.102. The mean metallic- 
ity measured with model A is in good agreement the value of 
Z = 0.27±0.15 determined by Edge & Stewart (1991), from 
EXOSAT observations of low-redshift clusters (using a sim- 
ilar spectral model. This value has been adjusted to account 
for the different Fe/H ratio assumed in that study.) For indi- 
vidual clusters, our results also show good agreement with 
those of Yamashita (1992), from GINGA observations, af- 
ter correction for the different definition of solar metallicity 
used in that work. The standaxd deviations in the metallicity 
measurements for the subsamples of CF and NCF clusters 
(5.7 X 10"^ and 6.5 x 10~^ for the CFs with models A and 
C, respectively, and 5.4 x 10"^ for the NCFs with model 
A) are smaller, by a factor of ~ 2, than that for the sam- 
ple as a whole, and are more in line with expectations from 
theoretical models {e.g. KaufFmann & Chariot 1997). 

We note that Abell 2319 is the only NCF-classified clus- 
ter with a metallicity comparable to that of a CF system. 
This is probably due to our conservative approach in classi- 
fying systems as NCFs. AbeU 2319 is the only NCF in the 
Table with a mean cooling time within 100 kpc of the cluster 
centre of less than 10^'^ yr. Further observations will clarify 
the identification of Abell 2319 as a CF or NCF cluster. 

Sixteen of the clusters included in our study (11 CFs 
and 5 NCFs) are also included in the study of Mushotzky 
& Lowenstein (1997). However, these authors only included 
a cooling-fiow component in the spectral analysis for two 
of the eleven CF clusters in their sample (Zwicky 3146 and 
Abell 2390). In general, our results show good agreement 
with those of Mushotzky & Lowenstein (1997), although our 
values for a few of the CF clusters are slightly higher. This 
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Figure 2. The metallicity measurements for the clusters (deter- 
mined using spectral model A) as a function of redshift. The figure 
illustrates a lack of evolution in the emission-weighted metallici- 
ties of clusters within z ~ 0.3. 



is presumably due to recent improvements in the calibration 
of the the ASCA instruments. 

We have searched for evidence of evolution in the clus- 
ter metallicities by splitting the clusters into two subsam- 
ples, with redshifts less and greater than the mean value 
of 2 = 0.210. This gives 16 clusters with z < 0.210 (mean 
redshift 0.14) and 14 with z > 0.210 (mean redshift 0.29), 
with an approximately oven mix of CF and NCF systems in 
the two samples (5 and 4 NCFs in the low and high redshift 
samples, respectively). The mean metallicities for the low— 2; 
and high— 2; samples (using spectral model A for all clusters) 
are 0.30 ±0.08 and 0.31 ±0.09, respectively. The application 
of the Students t-test shows the difference between the moan 
metallicities for the high— « and low— 2; samples to be signif- 
icant at only the 18 per cent confidence level. In agreement 
with Mushotzky and Lowenstein (1997) we therefore find no 
evidence for evolution in the emission- weighted metallicities 
over the range of redshifts covered by our study. We also 
note that we find no evidence for significant evolution in the 
X-ray gas mass fraction, with a mean value for both the low 
and high redshift subsamples of 0.16. 



4 DISCUSSION 

4.1 Metallicity gradients, mergers and mixing 

Fabian et al. (1994a) suggested that the higher emission- 
weighted metallicities for CF clusters could be due to inho- 
mogeneities in the ICM, with small blobs of cooler, denser 
and more metal-rich gas being immersed in hotter gas at the 
centres of the CF systems. The distributed mass deposition 
profiles within cooling flows requires that the ICM there is 
inhomogeneous (Fabian 1994 and references therein). Fabian 
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et al. (1994a) also noted the potential importance of abun- 
dance gradients in clusters. Abundance gradients appear to 
be a common feature of the cores of CF clusters. First no- 
ticed in ASCA (Fukazawa et al. 1994) and ROSAT PSPC 
(Allen & Fabian 1994) spectra of the Ccntaurus cluster, an 
increase in the iron abundance towards the centres of CF 
clusters has since been found in ASCA spectra for Abell 
496 (Hatsukade 1997), the Virgo cluster (Matsumoto et al. 
1996) and AWM7 (Ezawa et al. 1997). The abundance of 
iron across the Perseus cluster, which also hosts a large cool- 
ing flow, appears patchy with the highest value probably oc- 
curring in the core (Arnaud et al. 1994; see also Molendi et 
al. 1998). In contrast, ASCA studies of AbcU 1060 (Tamura 
et al. 1996), the Ophiuchus cluster (Matsuzawa et al. 1996) 
and the Coma Cluster (Watanabe et al. 1997), which have 
little or no cooling flows, show no abundance gradients. 

An abundance gradient will strongly influence the 
emission-weighted metallicity determined from the inte- 
grated X-ray spectrum of a cluster. The X-ray emission de- 
pends on the square of the ICM density and so is dominated 
by the innermost, densest regions of a cluster. This is partic- 
ularly the case for CF clusters, where a significant fraction 
(up to ~ 70 per cent) of the total X-ray luminosity may arise 
from within the central 2 — 3 hundred kpc (Peres et al. 1997; 
Allen et al. 1998). We have simulated the effects of metal- 
licity gradients on the mean emission-weighted metallicities 
for CF and NCF clusters using the observed X-ray surface 
brightness profiles for Abell 478 and 2218 as representative 
examples of CF and NCF systems. (These clusters have the 
highest-quality imaging data.) We find that a linear gradi- 
ent in metallicity dropping from ^ 0.8 solar at the centre to 
~ 0.1 solar at 500 kpc (and remaining roughly constant out- 
side this radius) leads to a mean emission-weighted metallic- 
ity, for a CF cluster like Abell 478, of about 0.4. For a NCF 
system, for which the X-ray surface brightness profile will 
be much loss sharply-peaked, the mean emission-weighted 
metallicity will be only ~ 0.2. Such gradients should be eas- 
ily detectable with AXAF. 

The abundance gradient model provides the most natu- 
ral explanation for the segregation in the metallicity results 
for the CF and NCF clusters. The metal-rich core may be the 
oldest part of a cluster, forming earliest in the deep poten- 
tial wells. This could enhance both the formation of massive 
stars and the retention of gas in these regions. (We note that 
many of the most-massive CF clusters also show evidence for 
ongoing star formation in their cores; e.g. Allen 1995). After 
their initial collapse, clusters continue to evolve by the ac- 
cretion of material, often via subcluster merger events. NCF 
clusters, such as those included in this study, are thought to 
have recently experienced a major merger event wherein a 
large mass component has strongly interacted with the clus- 
ter core (Allen 1998). Such events will significantly disrupt 
the X-ray gas in the core regions of the clusters and will mix 
and spread the central high-metallicity gas with the outer 
less metal-rich material. The abundance gradients in NCF 
clusters arc thereby reduced (or even destroyed), although 
the total mass of metals in the ICM is unchanged. (An im- 
plication of this is that the global metallicity of the ICM is 
more accurately estimated for the NCF systems, where the 
metals are more evenly mixed with the cluster gas.) The de- 
pendence of cluster metallicities on the presence or absence 
of a cooling flows then reflects whether these systems have 



had their core regions left undisturbed, or recently mixed, 
rather than on any internal property of the cooling flows 
such as their mass deposition rates. 

Reisenegger et al. (1996) have suggested that the abun- 
dance gradient in the Ccntaurus cluster could result from the 
cooling flow in that cluster concentrating the metals ejected 
by type-la supernovae in the outer regions of the cD galaxy. 
Although this mechanism is plausible for low-luminosity sys- 
tems like the Centaurus cluster, with relatively high ratios 
of the stellar mass in the central galaxy to the X-ray gas 
mass, it has more difficulty in accounting for the gradients 
in more X-ray luminous clusters, with lower stellar/X-ray 
gas mass ratios. This implies that the metallicity gradients 
in luminous CF clusters may have been present since some 
early epoch in their formation history. 

A further possibility to explain the metallicity gradi- 
ents in cooling-flow clusters is that most of the metals in 
cluster cores may reside in large grains. The lifetime of 
grains of radius ajj,m to sputtering in hot gas of density 
n is ~ 2 X lO'^a/nyr (Draine & Salpeter 1979). Provided 
that individual grains exceed 10//m in radius, they should 
survive for a Hubble time or longer throughout NCF ob- 
jects and beyond the cooling radius in clusters with cooling 
flows. Within cooling flows, the density rises inward so the 
grains are increasingly sputtered, releasing the metals into 
the gas phase which thus becomes increasingly metal rich 
toward the cluster centre. Such a model requires that about 
half the metals are originally injected into the ICM as large 
grains. We note that large grains are inferred to have formed 
and to carry most of the iron in the expanding remnant of 
SN1987A (Colgan et al. 1994). Within this model the typical 
mean metallicity in the core of a cluster would be 0.5-1.0 so- 
lar. The grains would also provide a possible source for the 
dust inferred in the central optical nebulosities in cooling 
flows (Fabian, Johnstone & Daines 1994b; Voit & Donahue 
1995; Allen et al. 1995) and may be related to the excess 
soft X-ray absorption observed in cooling flows {e.g. Allen 
& Fabian 1997). 

4.2 The effects of cluster evolution 

Within standard formation scenarios, clusters that form ear- 
lier are expected to have higher luminosities for a given tem- 
perature (since the gas density is higher at earlier formation 
epochs). Scharf & Mushotzky (1997) presented results from 
a single-temperature analysis of ~ 30 clusters observed with 
ASCA and demonstrated a positive correlation between the 
amplitude, Alt, of the Lx — Tx relation (where they deflne 
Lx = AltTx) and the mean (emission- weighted) metallicity. 
They suggest that the origin of this correlation is that clus- 
ters with larger values of Alt formed earlier and were better 
able to hold on to the metal-enriched gas expelled from their 
galaxies. Assuming that metals cannot be lost from the clus- 
ter potentials without a corresponding decrease in the X-ray 
gas mass, clusters with lower Alt should have lower baryon 
fractions by about a factor of two, given the observed range 
in metallicities. 

The Table lists the baryon fractions at a radius of 500 
kpc in the clusters, determined with the appropriate spec- 
tral models (spectral model C for the CF systems and model 
A for the NCFs). We see that the mean baryon fractions for 
the subsamples of CF and NCF clusters differ by only ~ 30 
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per cent, and thus that relatively little gas has escaped the 
potentials of the NCF clusters relative to the CF systems. 
(We note, however, that the distributions of baryon fraction 
values for the CF and NCF clusters are different. The appli- 
cation of a Student's t-test shows the mean values for the two 
subsamples to differ at > 99 per cent confidence. The NCF 
clusters also have a significantly smaller dispersion in baryon 
fractions than the CF systems, and a Kolmogorov-Smirnov 
test shows the two subsamples to be drawn from different 
populations at > 99 per cent significance.) If spectral model 
B rather than model A is used for the NCF systems {i.e. if 
the absorbing column density is included as a free parame- 
ter in the spectral analysis of the NCF systems) the mean 
baryon fraction for these clusters rises to 0.16, improving 
the agreement with the CF clusters. 

It is important to note that the clusters studied in this 
Letter are amongst the most X-ray luminous and by impli- 
cation most-massive clusters known. Material expelled from 
their member galaxies, even when part of an early low-mass 
subclump, is unlikely to escape from the total cluster po- 
tential. The accompanying Letter (Allen & Fabian 1998) 
discusses the impact of cooling flows on the LboI — Tx rela- 
tion for clusters. Allen (1998) discusses the effects of cooling 
flows on X-ray mass measurements. The new data presented 
in this Letter reveal a link between cooling flows and metal- 
licity measurements and suggest the presence of metallicity 
gradients in clusters with cooling flows. Such effects must be 
accounted for before attempting to determine cosmological 
parameters from X-ray observations of clusters. 
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